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SECTION  I 


INTRODUCTION 

The  general  approach  to  the  simulation  of  a nuclear  electromagnetic 
pulse  using  a radiating  system  consisting  of  planar  sources  has  recently 
been  considered  from  a qualitative  point  of  view and  methods  of  launch- 
ing the  required  transient  fields  are  currently  under  investigation.  Such 
planar  arrays  consist  of  aperture  sources  having  tangential  electric  fields 
which  are  excited  in  the  appropriate  time  sequence  and  spacial  distribution 
to  obtain  a desired  far-field  distribution.  Techniques  for  launching  the 
required  waves  for  the  simulation  of  an  electromagnetic  pulse  have  been 
discussed  in  previous  work  on  pulser  arrays! “ ’ 5 ' 

The  purpose  of  this  investigation  is  to  analyze  the  transient  radiated 
fields  of  planar  arrays  of  sources  which  might  be  used  to  simulate  a 
nuclear  electromagnetic  pulse  (EMP) . This  work  will  not  consider  the  problem 
associated  with  the  excitation  of  the  planar  sources,  and,  consequently,  the 
arrays  studied  will  consist  of  sources  with  known  aperture  field  distri- 
bution and  "time-delay"  distribution. 

The  first  phase  of  the  investigation  consists  of  the  analysis  of  the 
time-history  of  the  far-fields  of  an  array  on  N infinitely  long  slots  in  a 
perfectly  conducting  infinite  plane  as  shown  in  Figure  1.  The  analysis  of 
this  array  having  known  electric  field  distributions  is  based  on  a freon's 
function  approach  to  the  problem.  A two-dimensional  space-time  Green's 
function  for  a single  line  of  unit  amplitude  located  on  the  plane  is  pre- 
sented, and  the  actual  time  and  amplitude  distribution  for  each  source  is 
incorporated  in  the  solution  for  the  far-fields  through  the  use  of  appro- 
priate superposition  integrals  and  the  Green's  function.  The  related  in- 
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Figure  1.  Planar  array  of  two-dimensional  slot  sources. 
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vestigation  of  Sections  II  and  III  is  based  on  the  numerical  evaluation 
of  this  mathematical  formulation. 


SECTION  II 


TRANSIENT  RADIATED  FIELDS  FOR  A SINGLE  SLOT 
IN  AN  INFINITE  PLANE 


In  this  section,  the  transient  radiated  fields  of  a single  slot  are 
derived  as  a first  step  in  the  analysis  of  the  array  of  planar  slots  in 
an  infinite,  perfectly  conducting,  ground  plane  as  shown  in  Figure  1. 

I'hus,  the  idealized  problem  consists  of  a known  tangential  electric  field 
d istribut ion  in  the  plane  where  y=0  which  is  outlined  in  Figure  2.  For 
simplicity,  only  the  electric  field  component,  E^fx.t) , is  assumed  over 
the  slot  width  d,  and  the  tangential  electric  field  is  assumed  to  be 
zero  elsewhere  on  the  plane  as  a result  of  the  idealized  conducting  planes; 
however,  the  analysis  of  the  related  z component  or  the  combinations  of 
x and  z components  can  be  obtained  in  a similar  manner.  The  two-dimensional 
coordinate  system  used  in  the  mathematical  formulation  is  shown  in  Figure  3 
where  the  required  distance  factors  are  defined  for  the  assumed  source  and 
field  points.  The  medium  is  assumed  to  be  that  of  free  space,  i.e.  c = f(). 
u = u(1,  o = 0 for  y > 0. 

Because  the  tangential  field  distribution  is  known  everywhere  on  the 
y = 0 plane,  the  analysis  reduces  to  a problem  of  developing  mathematical 
expressions  for  the  direct  computation  of  the  transient  radiated  fields. 

One  approach  to  this  problem  is  to  invoke  the  surface -equivalence  theorem 
to  obtain  equivalent  surface  currents  from  which  the  desired  fields  can 
be  computed^.  ^ In  this  case,  from  Figures 2 and  3,  only  a magnetic  surface 
current  density  exists, which  is 
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Figure  2.  Geometry  for  slot  in  an  infinite,  perfectly  conducting 
ground  plane. 
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Figure  3.  Coordinate  Systems. 
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whore  F.  is  the  assumed  electric  field  distribution  existing  on  the 
in finite -length  slot  of  the  original  problem  and  the  e's  are  the  unit 
vectors  in  the  space.  The  equivalent  problem  for  y > 0 consists  of  an 
electric  conductor  in  the  plane  y = 0 (i.e.;zero  fields  for  y < 0)  with 
a magnetic  current  density,  J7.,  residing  on  the  surface  as  depicted  in 
Figure  4A.  This  equivalent  problem  can  be  formulated  in  terms  of  image 
theory  to  obtain  a solution  for  the  original  boundary  value  problem  where 
the  total  magnetic  current  density  is 


M = M 6(y)  = M 5(y)e  = 2E  e 

b L Z A Z 


Hence,  the  equivalent  problem  of  a sheet  or  ribbon  of  magnetic  current 
as  shown  in  Figure  4B  can  be  used  to  analyze  this  two-dimensional  problem. 

An  expression  for  the  transient  fields  of  the  strip  of  magnetic  current 
can  be  derived  in  terms  of  electric  vector  potential  from  Maxwell's  equa- 
tions for  time  dependent  fields.  Thus,  since  the  problem  involves  only  a 
source  of  magnetic  current  density,  i.evT=0,  p = 0,  p -0,  Maxwell's 
equations  for  this  problem  in  free  space  become 
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where  all  field  quantities  are  functions  of  both  time  and  space.  In  this 
homogeneous  media,  the  electric  field  can  be  represented  in  terms  of  an 
electric  vector  potential,  7, as 
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IT  = -v  x F 
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because  any  divergenceless  vector  is  the  curl  of  some  other  vector.  With 
this  vector  potential  representation,  the  magnetic  field  can  be  represented 
in  terms  of  F and  the  gradient  of  a scalar  as 


— R = ~ - V<j> 
£q  at  ym 


(8) 


where  <f>m  is  the  magnetic  scalar  potential.  From  the  previous  equations, 
it  can  be  shown  that  the  vector  potential  must  satisfy  the  wave  equation 


c2V2F  - ~ = -C2  M 5(y) 
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where  a Lorentz  gauge 
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1/2  _ 

is  assumed  for  c = 1/(vi0£q)  ' and  M = Ms(x,z,t)S  (y) . 

The  far  fields  can  be  computed  in  terms  of  the  electric  vector  po- 
tential from  the  previous  equations  if  equation  (9)  can  be  solved  in  terms 
of  the  known  magnetic  current  density.  In  this  case,  it  is  desirable  to 
reduce  the  potential  wave  equation  with  associated  boundary  conditions  to 
an  integral  expression  which  is  more  amenable  to  numerical  analysis.  Such 
an  integral  expression  can  be  obtained  through  the  application  of  the 
Green's  function  method  for  the  solution  of  differential  equation. 

The  problem  as  shown  in  Figure  4 reduces  to  a two-dimensional  problem 
in  free  space  where  the  solution  is  valid  for  y > 0,  and  equation  (9)  be- 
comes a scalar  wave  equation. 
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C2^l'z  - — r = -C2  Mz(x,Z,t)  6(y) 
3t 


because  from  equation  (2)  only  a z component  of  Ffs  exists.  Since  this 

is  a free  space  problem,  the  related  space-time  Green's  function  is  the 

principal  solution  of  the  general  equation, 


2 1 

Vx'v'  Gc  ' ~z  — = -5(x'-x»  y'-y»  t-t) 

1 c 3t 


cl 

which  has  the  solution 


(12) 


G0(x»y*t|x' ,y' ,T) 


f 0,  R>c(t-x) 

,R<c (t-t) 


c 

27 


(13) 


for 


(14) 


as  defined  for  a general  unit  line-source  in  both  space  and  time  in  Figure 
SA.  Special  boundary  conditions  at  infinity  and  causality  requirements 
dictate  this  solution  where  the  variables  of  (11)  have  been  interchanged  to 
insure  proper  function  dependence  in  the  integral  solution  of  equation  (11). 
This  Green's  function  is  the  characteristic  two-dimensional  solution  for 
the  scalar  wave  equation  caused  by  an  infinite  line  source  excited  at  r = t. 
Actually,  since  the  line  consists  of  many  point  sources  which  radiate 
spherically,  a response  at  some  distance  Rq  is  first  experienced  at  x = t 
and  subsequent  response  occuisfor  x's  > t;  thus,  the  response  decays  accord- 
ing to  equation  (13)  as  shown  in  Figure  5B. 
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Figure  5A.  Line  source  representation  for  a two-dimensional 
Green's  function  (pictured  for  y'=0) . 


Figure  5B.  Space-time  Green's  function  for  two-dimensional 
wave  equation  in  free  space. 
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Through  the  use  of  the  appropriate  Green's  theorem  and  initial 
conditions,  a solution  of  the  inhomogeneous  wave  equation  (11)  can  be 
exp i as  sod  in  terms  of  a superposition  integral  and  this  Green's  Function 
as 


fT  °°  r 

r~U.y,t)  = J j Mz(x' ,y' ,T)G0(x,y,t|x' ,y' ,T)dx'dy'dT 

0 -CO  — 


(15) 


where  the  required  frc*  space  boundary  conditions  and  causality  eliminate 
all  contour  contributions  and  T exceeds  any  time  t at  which  the  solution 
may  be  sought.  Now,  since  the  source  exists  only  in  the  x-z  plane  (y=0) , 
i.e.  M^(x,y,T)  = Mz(x,t)  6(y),  then  equation  (IS)  becomes 


f,(x,y,t)  = 


Mz(x'  ,-r)G0(x,y,t|x'  ,0,T)dx,dt 


(lb) 


0 -o 


which  is  an  integral  formula  for  computing  the  electric  vector  potential 
for  planar  two-dimensional  sources  of  known  electric  field  distribution 
through  the  use  of  equation  (2). 

The  time  dependent  or  transient  fields  radiated  in  the  far-field 
of  the  two-dimensional,  planar  source  can  be  expressed  in  terms  of  equa- 
tion (16)  with  the  aid  of  the  expressions  of  equations  (7)  and  (8).  Thus, 
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F = -V  x ez  j Mz(x’,t)  G()(x,y,t|x' ,0,T)dx'di 
0 -00 


(17) 


and 
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Mz(x\t)  G0(x,y,t|x'  ,0,x)dx'dT 
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(18) 


EZ_ 


- 


-r  v 


where  the  scalar  potential  component  of  the  H-field  equation  is  zero  in 
this  two-dimensional  problem.  In  the  three-dimensional  case,  this  term 
can  be  dropped  because  it  has  a higher  inverse  power  of  distance  than  the 
time  derivative  term,  and,  therefore,  the  scalar  potential  contribution 
is  smaller  than  the  vector  potential  contribution  to  the  field  at  large 
distances.  However,  the  scalar  potential  contribution  could  be  expressed 
in  terms  of  the  vector  potential  using  the  Lorentz  condition  of  equation 
(10)  and,  thereby,  included  in  a complete  three-dimensional  formulation. 

The  field  formulations  of  equations  (17)  and  (18)  are  not  readily  solved 
using  numerical  techniques  in  their  existing  forms,  but  for  some  simple 
magnetic  current  density  representations  these  expressions  can  be  reduced 
to  forms  with  tractable  numerical  solutions.  One  simplification  which 
greatly  reduces  the  numerical  computation  required  is  the  use  of  a step- 
function  dependence  in  time  such  that 


Mz(x,t)  = Mz(x)  u(t) 


(19) 


where  other  time  dependences  can  be  obtained  from  a superposition  integral 
based  on  the  summation  of  step  function  responses!^  For  this  case  of 
a planar  array  excited  at  t = 0 with  a step  function  dependence,  the  field 
expression  becomes 
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J GodT 
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dx' 
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(21) 


T 


where  the  functional  dependence  of  the  Green's  function  has  been  dropped 
for  simplicity.  In  the  present  form,  equations  (20)  and  (21)  re- 
quire double  integration  to  compute  the  field  quantities  which  is  usually 
undesirable  in  numerical  analysis;  however,  the  integral  in  the  dummy  time 
variable  of  the  Green's  function  can  be  obtained  in  a closed  form.  From 
equations  (20)  and  (13),  this  integral  becomes 


r = ^ 
G1  -V 


t-R/c 


0 


7c“(t-x)2-R2]1/2 


(22) 


and,  if  the  change  of  variables  n = t-x,  dn  = -dx,  and  R/c  = b are  used, 
then 


G1  = i- 


-dn 

TTTT77 


[ (n  -b  ) 


(23) 


Note  tl\at  the  upper  limit  of  (23)  is  evaluated  at  T = t-R/c  because  Gq=0 
for  T > t-R/c.  The  unit  step  response,  G^,  of  an  infinitely  long  line 
of  magnetic  charge  density  can  be  obtained  in  a closed  form  using  the 
indefinite  integral  form[^ 


dx/(x2-a2)1/2  = sin(x+(x2-a2)1/2) 

~1 


as 


G1(x,y,t|x')  = ^ 


in 


t + /t^~(R/c)‘ 



; t > R/c 


(24) 


Finally,  the  field  expressions  of  equations (20)  and  (21)  can  be 
written  in  terns  of  the  step  response  of  the  infinitely  long  line  (equation 
23)  as 

r°° 

(25) 


E = -V  x ez  Mz(x')  G^(x,y,t|x')dx' 
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and 
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Vl  = -em 


e.j  Mz(x')  G^(x,y,t |x')dx' 


(26) 


The  individual  field  components  can  now  be  obtained  in  cylindrical 
coordinates  defined  in  Figure  3 as 


r 


3G, 


E = M (x')  -^-dx' 
<j>  J z v J ap 

-oo 

3Gt 


E = - 
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and 


Hz  = 'e0 


3G, 

Mz(x')  ~dx’ 


(27) 


(28) 


(29) 


where  differentiation  under  the  integral  sign  is  valid  in  all  cases  be- 
cause the  limits  and  variable  of  integration  are  independent  of  the 
differentiation  variables.  If  the  chain  rule  for  differentiation  of  com- 
posite functions  is  used,  the  partial  derivatives  of  the  integrands  of  equa- 
tions (27)  and  (28)  can  be  evaluated  as 


3G1  8G1  SR 

K^x.y.t.x')  = — = ^ 


[t(t2-(R/c)2)1/2n2Up-x  cos^l 
2tt  R2(t(t2-(R/cV/2+t2  - (R/c)2] 


(30) 


and 


K2(x,y,t,x')  = 
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3R 
3 <p 


-[t(t2-(R/c)2)1/2+t2]x'sin» 

2ttR2 [t(t2-(R/c)2)  /2+t2- (R/c)2] 
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(31) 


where  G1  and  R are  differentiable  over  the  region  t > R/c  and  R and  p 
arc  defined  in  the  polar  coordinates  of  Figure  3 as 


R = [p2  sin2c|>  + (p  cosij>  -x')2]1^2 
p - [x2  . y2,1/2 


such  that 


3R  _ p-x'  cosiji 
3 P R 

3R  _ px'  sin<|i 
R 


Equations  (27),  (28),  and  (29)  can  now  be  expressed  in  simplified 
terms  as 


F:<t> " J Vx,>  ^1(x»y»t,x,)dxl 
- 00 

(36) 

ep  = ‘j  Mz(xr) K2Cx>y»t>x,)dx' 

- oo 

(37) 

and 

•OO 

Hz  = ' MzCx')  K3(x,y,t,x')dx' 

(38) 

where  K1  and  K2  are  defined  in  equations  (30)  and  (31)  and 

K3(x,y,t,x')  = aGl/3t  ' e0 
_ E0  1 

“ TT  — 7 TTTT  (39) 

(t  -(R/c)V/Z 

Equations  (36),  (37),  and  (38)  are  general  expressions  for  the  three 
far-field  components  produced  by  a two-dimensional  source  of  z-directed 
magnetic  current  density  for  an  assumed  step  function  excitation  in  time. 


; 


r-  ••>-*  • , ^ 


These  field  components  can  be  related  to  known  x-directed  electric  field 
distributions  through  the  use  of  equation  (2). 
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SECTION  I I 1 


SINGLE-  AND  MULTIPLE-  ELEMENT  ARRAYS  01'  SLOTS 
WITH  CONSTANT  SOURCE  DISTRIBUTION 

single  slots 

In  this  investigation,  slot  sources  of  known  constant  electric 
field  distribution  are  of  primary  interest.  Consequently,  as  a first 
step  in  the  solution  of  the  related  problem,  the  single  slot  with  an 
x-directed  electric  field  distribution  such  that 


?Ux')  = 2E0[u(x'+d/2)  - u(x'  -d/2)  ] 


(40) 


is  analyzed.  The  transient  field  components  for  the  single  slot  of 
figures  3 and  4B  with  constant  source  distribution  are,  from  equations  (36) 
through  (38)  and  (40) , 


(41) 


(42) 


and 


2,41 


h = 


0 z 


0/2 

2Zq  j K.(x,y,t,x’)dx’ 

-d  n 


(45) 


where  each  component  has  been  normalized  to  an  appropriate  constant  and  Z(1 
is  the  characteristic  impedance  of  free  space. 

The  expressions  of  equations  (41),  (42),  and  (43)  are  of  the  super- 
position or  convolution  type  integral  form  which  can  readily  be  represented 
in  graphical  form,  as  shown  in  Figure  6.  In  this  case,  these  problems 
consist  of  the  convolution  of  a constant  amplitude  pulse  function  in  space 
(x'l  and  the  appropriate  sharply  peaked,  kernel  function,  lntegrable 


r 
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Kernel  Function 


figure  6.  Graphical  representation  of  the  field  integrals 
for  a single  strip. 
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singularities  of  the  kernel  functions  do  appear  in  the  spatial  coordinate 
x'  for  t=R/c  where  R is  the  distance  between  the  field  point  and  source 
point.  Thus,  for  a given  point  in  space  (x,y)  and  a time  (t),  the  magni- 
tude of  the  field  is  related  to  the  total  area  underneath  the  curve  de- 
fined by  the  shaded  section  in  Figure  6.  For  the  general  case  (y  0),  a 
solution  which  increases  from  zero  to  a maximum  value  followed  by  a mono- 
tonically  decreasing  magnitude  as  a function  of  time  can  be  obtained  for 
the  eL  and  h field  components.  However,  the  expressions  for  the  e com- 

v Z p 

ponent  is  a direct  function  of  x'  and  sin  f>  (see  equation  31),  and  this 
component  vanishes  on  the  plane  y = 0 (J>  = 0,n).  In  addition,  this 
functional  dependence  can  change  the  polarity  of  the  integrand  of  equa- 
tion (42)  because  x'  is  an  odd  function  which  leads  to  field  expressions 

for  the  e component  that  are  either  positive  or  negative  depending  on 

P 

geometry  of  the  problem.  The  magnitude  of  the  time-history  of  the  e^  com- 
ponent is  similar  to  that  of  the  e,  and  h^  field  components  unless  the 
origin  of  the  coordinate  system  lies  within  the  slot.  The  single  slot  of 
Figure  4B  is  symmetrically  spaced  about  the  center  of  the  coordinate 
system,  and,  as  a result  of  the  ’’oddness"  of  the  integrand,  the  time- 
history  of  the  e component  has  both  positive  and  negative  excursions. 

The  rise-time  of  each  of  the  waveforms  of  these  field  components  is 
related  to  the  "delay-time"  between  field  contributions  propagating  from 
the  closest  and  furthest  source  points  on  the  slot  in  relation  to  the  field 
point  (See  Figure  3 for  geometrical  configuration).  This  rise-time*,  tf, 
can  be  defined  as  the  time  required  for  the  field  to  increase  from  zero  to 
one-hundred  percent  of  its  maximum  or  peak  value.  For  narrow  slots,  it  is 

a direct  function  of  the  width  for  e , e (in  all  but  the  special  cases), 

■f>  P 

*This  rise-time  definition  is  valid  only  for  single  slot  element  re- 
sponses which  increase  monotonically  as  a function  of  time  to  some  peak 
value . 
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and  h7  as  is  anticipated  from  the  "turn  on"  points  in  the  graphical  represen 
tation  of  spatial  convolution  in  Figure  6.  For  large  distances,  it  can 
be  shown  that  the  previously  defined  rise-time  for  these  field  components 
can  be  expressed  in  normalized  form  as 


Tr  = = (d/R(Pcos*  (44) 

where  R = R(| -d  cos  41,  tf  is  the  total  rise-time,  R(1  is  the  distance  from  the 
center  of  the  slot  to  the  field  point,  and  d is  the  total  width  of  the 
slot.  Equation  (44),  though  approximate,  clearly  indicates  that  the 
normalized  rise-time  is  a sensitive  function  of  the  angle  $,  particularly 
around  <f>  = n/2,  and  that  the  smallest  rise  times  are  obtained  for  broadside 
propagation  (<p  = n/2)  and  small  width- to-distance  ratios. 

.An  approximate  expression  for  late-time  behavior  of  the  et  component 
of  the  radiated  electromagnetic  at  large  distances  can  be  obtained  directly 
from  the  magnetic  field,  h . In  the  far-field,  i.e.  d/R  <<  1,  the  electric 

t*  u 

field  is  related  to  the  magnetic  field  by  the  characteristic  impedance 
of  free  space;  thus,  from  equations  (39)  and  (43), 


dx' 


, , 1/2 
[(ctr-fT] 


(45) 


when  t > R/c  such  that  a response  from  each  point  on  the  slot  is  experienced 
at  the  field  point.  Equation  (45)  can  now  be  expressed  as 


dx' 

, 1/2 

[a  + bx'  + cx'^1 


(46) 


which  is  intcgrablc  in  a closed  form,  where  from  the  geometry  of  the 
single  slot, 
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(47) 


2 2 

a = (ct)  - 


b = 2Rq  cosij) 


c = -1 


This  integral  has  a solution 
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a + bx'  + ex'  ] 


■ 1 -2cx-b 
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thus,  an  approximate  expression  for  the  electric  field  component  in  the 
t direction  is  obtained  in  the  form 


cos<J>+d/2R 

e sin  . 

4>  r / 2 i 

"Vcos-^+t^-I 


.j  cos<J>-d/2R0 
“^/cos^+t'-I 


where  = t/(R^/c),  and 


t 0 2 [ 1+ (d/RQ)  cosi))+  (d/RQ)  ] 


for  late-time  behavior.  Of  course,  equation  (49)  is  related  to  an  exact 
closed-form  expression  for  h based  on  equation  (45),  i.e.  e_  = Znh_;  and 
for  <j>  = it/2,  the  equation  (49)  reduces  to  a form 


2.-1 

e = - — sin  - 

<f>  " 7 

Vv1 


where  it  can  easily  be  seen  that  these  components  decav  in  time.  Note 

2 1/2 

that  for  [l+(d/lR())  ] in  equation  (51),  which  is  at  the  peak  of 
the  waveform, 


e.  = -1 


% 


for  the  special  case  when  $sir/2.  In  addition,  it  is  recognized  that 
equation  (49)  is  quite  general  and  can  be  employed  to  compute  the  complete 

late-time  behavior  after  the  peak  has  occurred  for  both  the  h and  cx 

z 4 

components  for  small  slot-width-to-distance  ratios. 

A numerical  solution  procedure  for  computing  the  time-history  of 
the  field  components  of  equations  (41),  (42),  and  (43)  has  been  developed 
and  programmed  for  a digital  computer.  This  procedure  is  a straightforward 
application  of  numerical  integration  using  an  open-ended  Gauss  quadrature 
formula  for  the  treatment  of  the  singularities  of  the  kernel  functions. 

The  application  of  numerical  integration  in  time-domain  problems  dis- 
cretizes the  representation  of  the  integral , i.e.,  the  total  field  is  the 
sum  of  weighted,  time-shifted  contributions  from  a finite  set  of  points 
on  the  slot  because  of  the  inherent  time-delay  between  integration  points. 
As  a result,  care  must  be  exercised  to  insure  that  the  spatial  variable 
of  integration  is  "sampled"  at  a sufficient  number  of  points  to  obtain  a 
convergent  integration,  especially  in  the  rise-time  region.  This  problem 
is  pronounced  for  small  observation  time  increments,  and  an  insufficient 
number  of  integration  points  results  in  an  irregular  waveform  consisting  of 
a serrated  shape  during  the  earlv-time  response  of  the  function. 

The  time-history  of  the  computed  field  components  using  this  numerical 
procedure  is  characterized  in  Figures  7 through  13  as  a function  of 
normalized  parameters 


T0  = (53) 

and  <j>/(n/2)  for  three  slot-width-to-distance  ratios,  d/R^, along  with 

approximate  fields  from  equation  (49).  As  can  be  seen,  each  e and  h 
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Figure  7B.  e vs.  Tn  for  d/R„  = 1/5  with  the  parameter  <J>/0/2) 


Figure  8A.  h_  vs.  t„  for  d/Rn  = 1/5  with  the  parameter  <t>/0/2)  > 0.75. 


Figure  8B.  h7  vs.  t,  for  d/R„  = 1/5  with  the  parameter  4>/0/2) 


510 


Figure  9.  e vs,  for  d/R-  = 1/5  with  $ as  a parameter. 


d/R  = l/30 


Figure  10B.  e vs.  tn  d/R0  = 1/30  with  the  parameter  $/(rr/2)  < 0.75. 


j>/[-rr/  2}  = 1.0 


Figure  11A.  h_  vs.  xn  for  d/R  = 1/30  with  the  parameter  */RQ  > 0.75 


d/R  = 1/30 
<*>/(tt/2)=0.75 


Figure  1-.  e vs.  tp  for  d/R(,  l/'>‘  with  f as  a parameter 


Approximate  Equation 


Figure  14B.  h_  vs.  t with  d/RQ  = 1/300  with  the  parameter  */(»/ 2) 


Figure  15.  e vs.  r.-,  for  d/R  , = 1/500  with  $ as  a parameter 


field  component  exhibits  a finite  rise-time  characteristic  followed  by  a 
decreasing  magnitude  as  a function  of  time  which  rapidly  approaches  a 
limiting  late-time  behavior.  Note  that  the  peak  of  the  functions  occurs 
later  in  time  and  tliat  the  peak  value  decreases  for  decreasing  values  of 
^/(fl/2).  This  "spreading  out"  of  the  waveforms  is  a result  of  the  time- 
delay  experienced  by  the  propagated  waves  from  the  slot  as  a function  of 
angle  which  is  related  to  the  "dispersion  effect"  in  the  frequency  domain. 

The  computed  transient  behavior  of  the  fields  of  equation  (41)  can  be  used 
to  illustrate  how  the  "spreading"  or  total  rise-time  of  the  waveform  varies 
as  a function  of  observation  angle  for  the  e^  component  of  the  single  slot 
as  shown  in  Figure  16.  Again  the  curve  emphasizes  that  the  rise-time  is 
dependent  on  angle  of  observation  and  that  it  changes  most  rapidly  near 
broadside  (4>=tt/2)  . 

It  is  interesting  to  note  that  the  numerical  results  of  Figures  7A, 

8A,  and  10A  indicate  that  the  time  response  of  a single  strip  has  a short,  but 
finite  rise-time  for  <t>/ (tt/20  = 1.0  (i  .e.j  broadside  to  the  slot).  This 
apparent  rise-time  is  a result  of  the  numerical  convolution  algorithm,  and 
it  can  be  shown  that  the  actual  time  response  has  zero  rise-time  in  this 

special  case  as  inferred  in  Figure  16.  From  Figure  16,  it  is  seen  that  the 

2 l /? 

early  time  response  for  <}>/ (tt/2)  is  determined  for  1<Tq< (1+  (d/2Rfi)  ) ; 

thus,  equation  (45)  reduces  to  the  form, 
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where  a = (ct)^  - Rg  and  0<a^^<  d/2.  This  integral  can  be  evaluated  as 


e.  = — 

d>  ir 


dp 
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(55) 
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d/Ro=l/30 

Approximate  rr,  Equation  (44) 


Measured  from  Graphs 
^ of  Time -History  of 
n e,  (d/R  =1/30) 


if  the  change  ol  variables,  u=x'/^a,  is  employed,  and  these  result  indicate 
tint  the  normalized  <j>  component  of  the  electric  field  is  independent  of 
time  during  this  interval  where  t'a’  is  finite,  however  small,  and  less 
than  d/2.  Therefore,  the  time  response  for  the  special  broadside  case 
has  zero  rise-time  and  a constant  magnitude  until  responses  from  all 
points  on  the  slot  have  reached  the  observation  point;  then,  the  field 
decays  in  a manner  described  by  equation  (51). 

It  can  be  seen  from  Figure  7B  that  the  approximate  expression  for 
the  <j>-component  of  electric  field  (equation  49)  is  in  error  for  4-/ (m/2) 

= 0.75  at  d/Rp  =1/5  while  the  h,  component  of  Figure  8B,  as  calculated 
with  the  approximate  equation,  is  exact  even  for  d/Rp  =1/5.  If  the  slot 
width-to-di stance  ratio  is  increased  to  1/30,  the  e^  components,  calculated 
by  both  numerical  and  approximate  expressions,  are  almost  identical  as 
shown  in  Figure  10B;  and,  if  d/R^  = 1/300,  it  is  difficult  to  distinguish 
between  the  two  computations  as  can  be  observed  in  Figure  13B.  Foi  '(-/2)=1, 
it  appears  that  the  approximate  equation  is  valid  even  at  d/R()  = 1/5  (See 
Figures 7A,  10A,  and  11A)  and  that  the  accuracy  of  the  approximation  is  a 
function  of  angular  position  as  w-ell  as  distance.  The  late-time  behavior  of  the 
h„  component  is  exact  in  all  cases  as  is  evident  from  the  graphs  shown 
in  Figures  8B,  1LA,  14A,  and  14B. 

Thus,  in  summary-,  it  has  been  shown  with  these  results  that  the  time- 
history  of  the  far-field  components  of  the  planar  slot  is  a critical  function 
of  both  angular  position  and  slot  width-to-distance  ratio.  The  c and 
h„  components,  which  have  the  same  basic  shape,  are  distorted  from  the 
peaked  broadside  response  for  angles  less  than  ninety  degree:  in  both 
amplitude  and  shape. 
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MULTI  ELEMENT  ARRAYS  OF  SLOTS 

With  the  results  already  obtained  for  the  single  slot  in  this  investi- 
gation it  is  relatively  easy  to  extend  the  theory,  as  well  as  the  numerical 
procedure,  to  the  array  problem  of  Figure  1.  Hence,  if  the  surface 
equivalence  theorem  is  applied  to  the  N-element  planar  array  of  slots  of 
Figure  17,  then  the  composite  magnetic  current  density  can  be  expressed  as 


M_(x,t)  = 2[En(x')[u[x'-(xn-dn/2)]-u[x’-(xn+dn/2)]]r(t- 


In)  (56> 


where  En(x')  is  the  electric  field  distribution  on  the  nth  slot,  x is 
the  distance  to  the  center  of  the  nth  slot,  t is  the  "turn-on"  time,  and 
d^  is  the  width  of  the  nth  slot.  Through  the  use  of  the  magnetic  current 
density  of  equation  (56)  and  the  theory  developed  in  Section  II,  equations 
(36) , (37) , and  (38)  can  be  used  to  obtain  integral  expressions  for  the 
iVelement  slot  array  of  the  form  (See  Figure  17  for  geometry) 
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vchere  causality  requires  that  (t-Tn)>R/c  for  each  individual  term  of  the 
series  expansion.  The  order  of  summation  and  integration  has  been  inter- 
changed in  the  formulation  and  the  kernel  expressions  are  the  same  functions 
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defined  in  equations  (30) , (31) , and  (39)  except  that  t is  now  replaced  by 
(t-x  ).  These  equations  still  represent  exact  expressions  for  the  fields 
of  this  two-dimensional  problem  involving  multiple  slots. 

Thus,  as  can  be  seen  from  these  expressions,  the  array  problem 
simplifies  to  the  summat ion  of  the  transient  responses  of  N single  slots 
with  the  appropriate  time  delay  and  position.  The  individual  transient 
response  for  each  of  these  slots  consists  of  a spat ial- shifted  convolution 
integral  of  the  type  investigated  for  the  single  slot  which  can  include  a 
time  delay,  t . A graphical  representation  of  the  convolution  integrals  for 
the  array  is  presented  in  Figure  18  for  two  times,  Tj  and  T?,  as  a function 
of  the  spatial  variable.  As  a result  of  the  spatial  distribution  of  the 
field  distribution,  it  is  easily  seen  that  the  sum  of  the  individual 
integrals  will  consist  of  time-delayed  responses  of  each  array  element. 

For  a finite  spacing  between  elements,  the  total  response  is  not  continu- 
ous as  a result  of  the  zero  contribution  between  array  elements;  and, 
consequently,  the  time-history  of  the  slot  array  is  characterized  by  a 
"serrated"  or  sawtooth-type  waveform.  The  shape  of  each  serration  of  the 
response  is  a direct  function  of  the  position  of  the  array  element  in  re- 
lation to  the  observation  point,  and  each  separate  contribution  can  be 
characterized  by  both  the  rise-time  and  delay  of  the  response  of  a single 
element  as  a function  of  $n/(ir/2).  The  time  delay,  t^,  between  peaks  of 
the  serrations  is  a function  of  the  spacing  between  edges  of  the  array 
elements  such  that,  from  Figure  17, 

Tsn  = W3  = CSn  ' (dnA-l)/2)cosVl/R0  (6°) 

where  is  the  normalized  time  delay,  Sn  is  the  center-to-center  spacing 
of  the  nth  ;ind  (n-1)^1  element.  For  an  array  of  equally  spaced,  constant 
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Figure  18.  Graphical  Representation  of  the  Component  field 

integrals  for  an  Array  of  Slot  Sources. 
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width  elements  the  time  delay  between  response  peaks  becomes 


Tg  = S cos<f>  (61) 

in  the  far-field  where  S is  the  normalized  element  spacing,  Sq/Rq,  and 
Sq  is  the  actual  element  spacing.  Thus,  the  slot  array  is  completely 
characterized  by  the  response  of  the  individual  array  element , the  spacing 
between  elements  and  the  total  number  of  the  element  which  determines  the 
aperture  size  and,  therefore,  the  region  of  maximum  response. 

The  expression  for  the  ^-component  of  electric  field  of  the  N-element 
array  (equation  57)  has  been  solved  using  the  numerical  procedure  utilized 
for  the  investigation  of  a single  slot  source.  The  results  of  this  computer 
analysis  for  a two  and  a three-element  array  of  equally  spaced,  constant- 
width  elements  are  presented  in  Figures  19  and  20,  respectively,  for 
different  element  spacings.  For  S=1.0,  i.e.,  all  elements  touching,  the 
responses,  with  symmetrical  placement  of  the  array,  are  simply  the  transi- 
ent response  of  a single  slot  element  of  either  twice  or  three  times  the 
element  source  width  for  the  appropriate  array.  This  response  consists  of 
a unit  magnitude,  single  pulse  with  zero  rise-time  and  a long  "fall-time” 
characteristic  of  the  single  slot  responses  of  Figures  7,  8,  and  10  for 

. For  S>1,  the  time -history  of  the  total  array  consists  of  the  sum 
of  the  time -delayed  responses  of  each  of  the  array  elements.  For  the  three - 
element  array,  the  response  separates  from  the  S=1  case  into  a waveform  having 
two  peaks.  One  of  these  peaks  is  a response  from  the  center  slot  and  the 
other  results  from  the  remaining  two  symmetrically  spaced  elements  which 
arrive  at  the  same  time.  There  is  a striking  difference  between  these 
contributions  because  the  response  from  the  center  element  is  not  drastically 
distorted  because  the  time  delay  across  this  element  is  small  and  the  time- 
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delay  for  points  on  the  off-center  elements  increases  as  a function  of 
position  which  results  in  deformation  of  the  basic  shape  of  the  response. 

As  shown  in  the  other  curves  of  Figure  19,  the  second  peak  is  seen  to  decrease 
in  both  magnitude  and  width  and  to  occur  at  later  times  for  increases  in 
element  spacing.  Numerical  results  for  the  two-element  array  are  similar 
to  that  of  the  three  element  array  without  the  first  peak  of  the  center 
element,  as  would  be  expected.  The  transient  responses  of  5,  7,  and  9- 
element  arrays  of  equally  spaced,  constant -width,  slot  sources  is  shown 
in  Figures  21,  22,  and  23,  respectively,  for  (|>(it/2)  equal  to  1 which 
shows  that  the  array  response  is  indeed  the  sum  of  the  time-delayed  re- 
sponse of  each  array  element.  In  addition,  the  effects  of  changes  in 
element  spacing  and  element  widths  in  a 9-element  array  are  dramatized  in 
Figures  24  and  25. 
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SECTION  IV 


CONCLUSIONS 

This  report  presents  general  expressions  for  the  transient  radiated 
electric  and  magnetic  fields  for  a planar  array  of  sources  with  elements 
having  known  electric  field  distribution,  as  well  as  time  delay,  located 
on  an  infinite,  perfectly  conducting  plane.  The  analysis  of  such  arrays 
is  based  on  the  space-time  Green's  function  approach  to  the  electric  vector 
potential  problem  rather  than  the  more  common  frequency  domain  approach 
applied  to  related  problems! ^ 

Planar  arrays  of  two-dimensional  sources  are  analyzed  using  numerical 
computation  of  the  appropriate  superposition  integrals  in  conjunction  uith 
related  Green's  function  for  constant  amplitude  distribution  over  the  in- 
dividual apertures.  Plots  of  the  time  history  of  the  radiated  far-fields, 
as  well  as  near  fields,  are  presented  for  a single  element  and  other  array 
configurations  for  the  constant  amplitude  distribution.  Results  of  the  fielc 
amplitude  variation  as  a function  of  observation  angle  are  presented  for  the 
single  element.  Note  that  the  time  history  of  an  element  is  extremely  sensitive 
to  observation  angle  because  of  the  time  delay  induced  by  the  change  in  distance 
a given  response  must  travel  to  reach  the  observation  point. 

The  time  history  of  the  fields  for  these  array  elements  is  character ist 
of  the  two-dimensional  wave  equation  solution  where  the  history  consists 
of  the  response  from  the  point  on  the  element  nearest  the  point  of  observation 
followed  by  responses  which  arrive  later  in  time  from  points  more  distant  on 
the  infinitely  long  aperture.  It  is  evident  from  the  responses  of  the  multi- 
element arrays  tliat  the  total  response  of  the  system  is  simply  the  I incar 
superposition  of  the  response  of  each  element  appropriately  delayed  in  time  for 
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the  existing  observation  angle.  in  addition,  from  Figure  24,  it  can  be 
observed  that  the  time  response  at  large  distances  from  a finite  array 
approaches  the  response  of  a single  element  for  broadside  observation  because 
time  delay  to  each  element  approaches  zero. 
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